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04 ' Glauber theory for nucleus-nucleus scattering at high incident energies is reformulated 

f-^ ' so as to also be applicable to the scattering at intermediate energies. We test the validity 

ij N , of the eikonal and adiabatic approximations used in the formulation, and discuss the rela- 

ry\ . tionship between the present theory and the conventional Glauber calculations using either 

the empirical nucleon-nucleon profile function or the modified one including the in-medium 

^ ' effect. 

(N 

§1. Introduction 






Experiments with radioactive beams of unstable nuclei have opened new frontiers 
in nuclear physics. New features of unstable nuclei such as a halo structure were 
P^ I revealed; see, for example, Ref. 1). The Glauber theory^) has widely been used as a 

powerful tool for studies of reactions with unstable nuclei observed at intermediate 
ff^ . energies such as 50 - 800 MeV/nucleon.'^'''^) It was reported,^-*' ^' however, that 

^ I some modifications of the nucleon-nucleon (NN) scattering profile functions in the 

eikonal approximation were necessary in order to reproduce the data at energies 
[~^ . of less than 500 MeV/nucleon.^^'^' Such phenomenological modifications obviously 

f^ I require theoretical foundations. In this paper, we address this problem. 

t^ ' The Glauber theory describes the scattering of two nuclei P and A as collisions of 

all nucleons in P with those in A. The theory starts with the many-body Schrodinger 
^ ] equation 

[E-K-hp-hA-V]'F = 0, (1-1) 

'k>( I where V = YlieP 7ga%' '^ith Vij the NN interaction potential, E is the energy of 

^ ■ the total system, K is the kinetic energy operator of relative motion between P 

and A, and /ip (/ia) is the internal Hamiltonian of P (A). Assuming the adiabatic 
approximation for the internal motion of P and A and the eikonal approximation, 
the theory gives the scattering amplitude at high energies and small scattering angles 
as 

//3a = ^ / db e*-?" (^^|1 - n(l - rNN(b.,))|<^a), (1-2) 

i,j 

where hk (hq) is the initial (transferred) momentum, b is the component of the 
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relative coordinate R between the centers of mass of P and A perpendicular to fc, 
^a (^/s) is the antisymmetrized internal wave function of the initial (final) channel, 
and rij = {bij, Zij) is the displacement of i from j with 6,-,- (zij) the component of Vij 
perpendicular (parallel) to k. Here, /NN(&ij) is the profile function of the scattering 
of nucleon i in P and nucleon j in A given by 

Av[N(feij) = 1 - exp / dzij Vij{rij) , (1-3) 

where v^-ei is the relative velocity. 

If the eikonal approximation is valid for NN scattering in free space, r^-^{hij) in 
(jl-3p should agree with the Fourier transform r^ of the NN scattering amplitude 
/nn determined from the data on NN scattering in free space, 

r^'^ih,) = ^ /e-9^.-^»^/NN(Q,,)dg,,-, (1-4) 

where kij (Qjj) is the initial (transferred) momentum of relative motion between 
the two nucleons i and j in free space and the 2-dimensional integration is over the 
components of q^j perpendicular to kij. Instead of using InnC^ij); -HInI^u) ^^^ 
been used customarily, since first introduced in Ref. 7). However, it should be noted 
that the replacement of Inn(^«j) by r^'^{bij) in (|l-2p is correct only if the eikonal 
approximation is valid for the NN scattering in free space, since (|l-3p has already 
been used in the derivation of (jl-2p . 

Substituting (jl-4p into ()l-2p . one can calculate the nucleus-nucleus (PA) scatter- 
ing amplitude //Jq. In Ref. 7), /nn(Qjj) was given, neglecting the Coulomb effects, 
in the form 

/NN(g.,) = /NN(0)e-^'9-/2 (-^.5) 

with 

fNN{0) = {i + a)kija/4, (1-6) 

where a is the NN total cross section, and the constants a and /3 were derived from 
the experimental data.^' This prescription works well at high energies. As already 
mentioned, however, (jl-6p had to be modified at energies lower than 500 MeV. The 
parameters a and (3 had to be adjusted away from those made to fit the PA scattering 
data. This cast doubt on the prescription, in particular on the validity of the eikonal 
approximation to NN scattering. In fact, the condition for the validity of the eikonal 
approximation to NN scattering is 

\vij/Eij\ < 1, kijtt > 1 , (1-7) 

where Eij and kij are, respectively, the kinetic energy and the wave number of the 
relative motion of i and j, and a is the width of the region(s) in which Vij changes 
rapidly. However, since Vij has a strong short-range repulsive core, for example 
Vij ~ 2000 MeV at rij = in the case of the realistic NN potential AV18,^) it is 
obvious that the first condition of (jl-7p is not satisfied. 

In order, therefore, to examine the accuracy of the eikonal approximation to NN 
scattering, we calculate the NN scattering amplitude at the laboratory energy E'nn = 
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300 MeV with the eikonal approximation and compare the result with the exact one. 
Figure 1 shows the on-shell NN scattering amphtude /NN(Qii)- For simphcity, we 
take only the central part of the realistic NN potential AV18 for the triplet-even 
state. The solid (dashed) and dotted (dash-dotted) lines show, respectively, the 
real and imaginary parts of the resulting scattering amplitude of the exact (eikonal) 
calculation. The eikonal amplitude deviates considerably from the exact one, even at 
small q, for both the real and imaginary parts. In particular, the deviation is serious 
on the imaginary part. This means that the Glauber theory can not accurately 
predict the reaction cross sections of nucleon-nucleus (NA) and nucleus-nucleus (PA) 
scattering. This test clearly shows that the use of p-4p in (|l-2p is inaccurate. 
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Fig. 1. Test of eikonal approximation for on-shell NN scattering amplitude /NN(q) at the laboratory 
energy -Enn = 300 MeV. As the NN interaction, the central part of AV18 for the triplet-even 
state is taken. 



A systematic way of making non-eikonal corrections was proposed in Ref. 9) for 
high-energy potential scattering. In the method, the transition matrix is expanded 
in a power series oi S = \vo\ / (khvrei) , where vq is a typical potential strength. The 
method, however, has not been applied to intermediate-energy NN scattering because 
6 > 1 in such a case. 

For the NA scattering mainly at high incident energies, some methods for treat- 
ing non-eikonal effects have been proposed, ^'^''^^' but they require much more diffi- 
cult and/or complicated calculations than the Glauber approximation. For example, 
Wong and Young introduced the pseudopotential f^ • that reproduces the empirical 
NN scattering amplitude in the eikonal approximation and estimated the corrections 
to the Glauber amplitude when u^- is used instead of Vij}^' However, it is not easy 
to evaluate the corrections induced by the difference Vij —v'-, since \v[-\ <C |%| and, 
therefore, Vij — v[; is much larger than v[,. 

It was proposed in Ref. 5) that f-^]<f{k',k) is modified from the empirical NN 
scattering amplitude in free space to that calculated with the Brueckner (7-matrix. 
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Although this seems a reasonable proposition, it is obvious that the theoretical foun- 
dations of such a phenomenological procedure need be examined. 

In this paper, we propose an accurate and practical method of treating the 
intermediate-energy PA scattering within the framework of the Glauber theory. For 
this purpose, we use the effective NN interaction r that describes the NN collisions 
in PA scattering instead of the bare potential Vij. We introduce this effective NN in- 
teraction in the Glauber theory on the basis of the multiple scattering theory (MST) 
of Watson^^^ using the formalism of Kerman, McManus, and Theler (KMT).^^) It 
turns out that the r satisfies the condition (|l-7p better than Vij and can even take 
account of the modification of the NN interaction in the nuclear medium that has 
long been thought to be necessary for reactions at low and intermediate incident 
energies. ^^ 

For the high-energy nucleon-nucleus scattering in which the Glauber approxima- 
tion is good and, therefore, ()l-3p is accurate, the relationship between the Glauber 
theory and MST was investigated in detail. ^^^'^^^ These studies, however, were 
mostly concerned with the cancellation in the Watson expansion between the re- 
flection terms and the off-pole contribution of nonreflective terms. ^^' This is not 
addressed in the present paper. We present a method of going beyond the ordinary 
Glauber theory by calculating all the Watson series using the Glauber approxima- 
tion. 

This paper is organized as follows. In §[2l we present a new version of the Glauber 
theory based on the effective NN interaction r. In § [3l some important points are 
discussed. In § 13.11 we approximate r as a two-body operator such as the transition 
matrix of NN scattering in free space or the Brueckner ^(-matrix, and we presents a 
way of localizing the two-body operator. In § 13.21 we test the eikonal and adiabatic 
approximations used in the present formalism. Such a test is not feasible for PA 
scattering consisting of an infinite series of multiple NN scatterings. The test is then 
made for a single NN collision in PA scattering, which is an essential building block 
of PA scattering. In ^ 13.31 we discuss the relationship between the present theory and 
the conventional Glauber calculations using either the empirical NN profile function 
(|1-4|) or the modified one including the in-medium effects.^' Section H] is devoted to 
a summary. 

§2. Formulation 

2.1. Multiple scattering theory for nucleus-nucleus scattering 

The original KMT formalism-'^^' is for nucleon-nucleus (NA) scattering. We first 
extend it to the case of nucleus-nucleus (PA) scattering between P and A. The 
transition matrix of PA scattering is given hy T = V{1 + GqT) with 

V 

^0 = E-K-h^-h^^ie' ^^'^^ 

where V = PpPa with V-p {'Pa) the projection operator onto the space of antisym- 
metrized wave functions of P (A). 
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The transition matrix T is given by 

T= Y. ^*^' (2-2) 

iePjeA 

where Tij satisfy a set of coupled equations 

Tij = Tij + Tij&o< y ^ Tmn ~ Tij > , (2"3) 

meP,neA 

where 

Tij=Vij{l + GoTij) (2-4) 

is an operator that describes the scattering of nucleons i in P and j in A. The 
potential Vij, in general, contains operators acting on the spins and isospins of i and 
j, which we suppress for the simplicity of notation. 

The proof of (|2-3p is as follows. We define X as X = J2ieP /eA-^u- Equation 
([2^ is then reduced to Tij = {l + TijGo)~'^Tij{l + GoX). Identifying {l + TijGo)~^Tij 
with Vij in the equation and summing Tij over i and j, one can obtain X = V{1 + 
GqX). We then find that X = T. 

Because of the antisymmetry of the nuclear wave functions, which is maintained 
by V, the matrix elements of the operators Tij and Tij are independent of the labels 
i and j, so that (|2-3p can be written as 



= T + {Y-l)TGoe, (2-5) 

where 6 = Tij and Y = Np x A'a, with A'^p (Aa) the nucleon number of P (A). 
Multiplying (j2-5p by y — 1, one obtains 



T' = U{1 + GqT') , (2-6) 

where 

T' = {¥ - 1)9, U={Y -1)t . (2-7) 

Then, T is obtained from T' as 

Y 

T = T' . (2-8) 

Y -1 ^ ^ 

Equations (j2-6p and (j2-8p constitute the final results of the extended KMT. The 
antisymmetrization between the incident nucleons in P and the target nucleons in A 
has been neglected so far. It was shown, ^^^'^^^ however, that it can be taken care of 
using r which is symmetrical with respect to the exchange of the colliding nucleons. 

2.2. Glauber approximation 

We proceed to the calculation of the matrix elements T'^^ for the transition 
a —^ P using the Glauber approximation. We first define the wave matrix f2a that 
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gives the wave function with an incident wave in channel a by iP'q = Qa ^a4>a, 
where (pa = (27r)~^''^ exp[ifc • R]. The wave function ^a satisfies 

{K + hp + hA + U- E)'^W = . (2-9) 

The matrix elements of T' for the transition q — > /? are then given by TL^ = 

{^l3(j)p\U\'I'a )■ In the present formalism, (Pa is an eigenstate of the internal Hamil- 
tonian hp + h\ with the realistic NN interaction Vij . In actual calculations, however, 
the Vij is often replaced by an effective interaction that usually has a weak repulsion. 
The short-range correlation neglected by the replacement may not be important at 
forward PA scattering with small q. 

The potential U can be written in the form 

C/=^ E ^^^ (2-10) 

ieP, jeA 

because of the total antisymmetry of the wave functions of P and A. By definition 
(j2-4p . Tij is a many-body operator acting on all the nucleons of the total system. 
We assume, however, as commonly performed in practical applications of MST and 
discussed in detail in S I3.ll that Tij can be well approximated by a two-body operator 
acting only on i and j and in coordinate representation by a local potential depending 
only on Vij. 

Equation (|2-9p then has the same form as the original Schrodinger equation 
(jl-ip for the wave function of the total system except that V is replaced by U. The 
Glauber approximation can, therefore, be applied to (j2-9p if the same conditions are 
satisfied, i.e., the adiabatic approximation to the internal motion of P and A is used 
to approximate /ip and /ia by their ground-state energies, and the conditions for the 
eikonal approximation, 

\{Y -1)t/E\<^1, A;a^»l, (2-11) 

are satisfied, where a,- is the range of the region in which r changes rapidly. Tij has 
a much weaker dependence on rij than Vij, as shown in S I3.1i Hence, it satisfies the 
conditions in (j2Tip better than the Vij of the original Glauber theory. 

We denote the center of mass of P (A) by Rp (Ra) and the coordinate of 
nucleon i (j) in P (A) by rj (rj). Then, R = Rp — Ra is the relative coordinate of 
P from A by which the scattering is described, and Xi = Vi — Rp {yj = Vj — Ra) 
is the intrinsic coordinate of i (j) in P (A). With this notation, U = U{R,(,) where 
^ = ({ajjjp, {y Ia). For simplicity of notation, however, the intrinsic coordinate ^ 
is suppressed in the following. 

In the Glauber approximation, TL are given by 



T^, = C j dhe"i-''{'Pp\ru{h)\<Pa) , (2-12) 



where C = —ifi? k / {{2tt)'^ jia) with /u^ the reduced mass in the initial channel a, and 
the profile function of PA scattering is given by 

rf/(6) = l-exp[ixt/(6)], (2-13) 
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Xu{b) = -^ r dzU{z,h), (2-14) 

where v^ei = hk/fia- Using (j2-10p . one can rewrite the phase shift function xu{b) as 

Xu{b) = ^ ^ xS?(bi,), (2-15) 
iGP, jeA 

where 

Xnn (^y) = --^^ / dzij T{zij,bij) (2-16) 



oo 



is the phase shift function of NN scattering by effective interaction r. The transition 
matrix elements are then given by 

Y 



'-(3a 



-C I dhe'i-''{^p\ru{b)\^^), (2-17) 



Y 

P A 



Fuih) = 1 - n 11(1 - r^^\h,)), (2-18) 

r^^\h,,) = 1 - exp[^ixS? (&,,)]• (2-19) 

The scattering amphtude fpa and the cross section a^a are given by 

respectively. Equations (j2-17p - (j2-19p with (j2-16p constitute the principal result of 
this paper, where the conventional form of the Glauber theory with (|l-2p and p-4p 
is reformulated using the effective interaction r. 

Most PA colhsions satisfy F > 1. The factors {Y - l)/Y and Y/{Y - 1) in T 
only have appreciable effects on collisions between very light nuclei. For the other 
collisions, the resulting scattering amplitude has the form of the original Glauber 
theory except that the bare NN interaction v is replaced by r. 

2.3. Elastic scattering 

For elastic scattering, where a = (5, one can rewrite (j2-17p in the form 

r„„ = TF^C / dh e'''-^(l - e^>^°p'('')) (2-21) 

with 

Xopt(&) = -iln(^„|exp[ixc/(6)]l'?'a) . (2-22) 

Since Xopt(^) is a function of b, there must be a one-body potential Uopt{z-,h) such 
that 

1 /•°° 
Xopt(b) = -J—J dz C/opt(2, b) . (2-23) 
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Since Uopt{z,h) is a one-body potential that describes the elastic scattering, it is 
the optical potential. This definition of Uopt{z,h) differs from the ordinary one by 
the overall factor Y/{Y — 1) in T, which is, as has already been mentioned, nearly 
equal to 1 except for scattering between very light nuclei. The potential Uopt{z-,b) 
is obtained as the solution of (|2-23p : 

One can obtain T^o by solving the Schrodinger equation using the optical potential 
f/opt(r-). 

In the case of |xc/(^)| < 1) exp[ixc/(&)] in (|2-22p can be expanded in powers of 
Xu{b). This leads to the cumulant expansion of Xopt(&)-^ If one retains only the 
lowest-order term, one obtains the optical limit, Xopt(^) ~ {^a\xu{b)\^a), and 

<i(^) = ^(^«l E r{z.,,Mi)\<l>a) , (2-25) 

ieP, jgA 

where the superscript (0) stands for the optical limit. Equation (j2-24p includes all 
orders of the cumulant expansion, although the calculation is not easy. 

The optical theorem citot = 4:TTlioafaa{k,k)/k for the total cross section Utot 
yields from (|2^ and (|2^ 

2Y f 

cTtot = y—^ I db Re{'Pc,\ru{b)\<P^) . (2-26) 

The angle-integrated cross section of elastic scattering is 

(Tel= / df2k' \faa{k',k)\^ 

^ {^y 1^^ \{^^\ru{b)\^a)\' . (2-27) 

When the elastic scattering is concentrated in the forward direction, df2k/ is nearly 
on a plane perpendicular to the direction of k, so that df2j^i ~ d^k' /k' , where d'^k' 
is an area element on that plane. ■^^ The reaction cross section Crcac is then obtained 
as 

2 



f^reac = C^tot " ^el = dh [l - \Z\' , (2-28) 



where Z = -1/{Y - 1) + Y/(Y - 1) • S'^^ with 

5^„ = (<?<, |exp[^ E iX^NN ibij)]\'Pa). (2-29) 

i€P,jeA 

In the case of 1" ^ 1, which most collisions satisfy, the reaction cross section of 
(|2-28p is reduced to that in the ordinary Glauber method with v replaced by r. 
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§3. Discussion 

3.1. Effective interaction 

The effective interaction r is the key to not only the present formalism but 
also to various theories of direct reactions. The effective interaction r is a many- 
body operator, because Eq. (j2-4p for r includes the many-body Green's function 
Go, which depends on all the internal coordinates of P and A. There is a long 
history of development of various approximations for r. Those widely used in many 
practical applications are to simply replace r by a two-body operator. In the impulse 
approximation, r is approximated by t, the transition matrix of NN scattering in 
free space. ^'^^ This has been successful for describing NA scattering at high energies 
of E^A > 500 MeV.i^)'^*^)'^^) At lower energies of E^a < 500 MeV, the effects of the 
nuclear medium on r become significant,^^'' ^^)'^^' and the Brueckner (^-matrix has 
been used for NA and PA scattering. ^''^"^^^ 

The t and g matrices are nonlocal operators. In many of their applications to 
the analysis of experimental data, however, they have been given as local potentials 
in coordinate representation. Love and Franey^ ' presented the t matrix elements in 
the form of superposition of Yukawa potentials that reproduce on-shell elements of 
the t matrix. In many of the g matrix applications, it is calculated for nuclear matter, 
parameterized in the form of the superposition of Yukawa or Gaussian potentials that 
reproduces on-shell or half off-shell elements of the g matrix, and translated to finite 
nuclei by a local density (LD) approximation. Although the justification of the LD 
approximation remains as a fundamental question, the use of g as r is successful in 
reproducing the NA scattering over a wide energy range of 50 < -Ena ^ 800 MeV.^^^ 
For PA scattering, however, there is an ambiguity in the definition of LD.^'^' Although 
further progress is required on this point, the g matrix may be one of the most 
plausible substitutes for r. 

The local potentials thus obtained are ambiguous in their rij dependence partic- 
ularly at small rij. Actually, the Yukawa- type local potential is singular at rij = 0, 
but the Gaussian-type is not. Obviously, the singularity is an artifact. Since the 
accuracy of the eikonal approximation is sensitive to the dependence, the localized 
potentials mentioned above cannot be used in the present theory. 

In order to determine a local r matrix from a nonlocal one T{rij,r^-) with no 
ambiguity, we introduce the trivially equivalent local (TEL) r matrix t (vij) as 

v{r,J)^P = Ti°^(ri,)0o (3-1) 

with (j)Q = {27r)~^''^ ex.p{ikij ■ rij). In the case of r = t, ■0 is the solution of the 
Schrodinger equation of the NN system with v, while in the case of r = 5 it is the 
solution of the Bethe-Goldstone equation. Since the discussion is parallel between 
the two cases, we consider, hereafter, the case of r = t. 

The TEL r matrix, t (rij), is related to the nonlocal matrix as 

r'°%r,,)Mr^J) = [ dr^r(r,„ r^)0o(r-:,) , (3-2) 
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because the left-hand side of (13-ip is equivalent to the right hand side of (l3-2p . The 
Fourier transform of the right-hand side of (j3-2p gives the half off-shell elements of 



T{rij,r'--). Equation (j3-2p ensures that t (rij) gives the same half off-shell elements 



U ^ ' ij) 



as the nonlocal matrix T{rij, r'-). Thus, T'°'^(rjj) is derived from the nonlocal matrix 
with no ambiguity. 

Even iiv{rij) is a central force, r'°'^(rjj) is not a central one, because it depends 
on the angle 9 between Vij and the initial momentum hkij (the z axis). This is not 
a problem in the formalism of § EJ since the integration of r 



loc. 



Ti 



^loC^L \ _ I J„. _l0C 



f'°%bij) = J dzi,T'°'irij), (3-3) 

is an input. The 9 dependence of T^°'^{rij) is integrated out in f'°'^(6jj) through the 
integration over Zij = rij cos{9). 

Usually, the TEL potential r is a function of the relative coordinate Vij and 
Vr-j^ , and the spin and isospin of the colliding pair. When the target is a double- 
magic nuclus, the tensor part of r does not contribute to U^ [. For this reason, in 
most analyses based on the folding model, only the central and spin-orbit parts are 
taken into account. Actually, the tensor force is difficult to include in the Glauber 
theory because of its noncommutative character.^) Thus, we assume 

ri°^ = T, + TsoS ■ L , (3-4) 

where L = —ivij x Vr-^ and U is the sum of the Pauli spin matrices for the two 
nucleons i and j. The angular momentum hL can be replaced by hvij x kij if 
\Tso/Eij\ <C 1 and kijago ^ 1, where aso is the width of the region in which Tgo 
changes rapidly. The potentials of the form (j3-4p are most suited to the Glauber 
approximation, which we assumed in § [2l 

Now we calculate f by solving the Schrodinger equation of the NN system 
with V. For simplicity, we neglect the spin and isospin of the colliding pair. As 
v{rij), we take the central part of AVIS^-* for the triplet-even state. The results 
are shown in Fig. El which represents the Zij integrated TEL potential f as a 
function of bij for the laboratory energies of £'nn = 300 and 800 MeV. The solid 
(dotted) curve is the real (imaginary) part, while the dashed curve corresponds to 
v{bij), the Zij integration of v{rij); at bij = 0, v = 1587 MeV ■ fm and f^°'^ = 
278 - 86z (267 - 212i) MeV • fm for ^nn = 300 (800) MeV. As expected, f i°'= has 
a much weaker bij dependence than v. The TEL potential r thus obtained is not 
singular at bij = 0. It is natural to think that the TEL potential of the g matrix also 
maintains the same property. 

3.2. Validity of eikonal and adiabatic approximations in derivation of F^^ 

The effective NN profile function F^-^ (bij) defined by (j2-19p . which describes 
a collision of the i-j pair in the PA scattering, is the key element of the present 
formalism, as shown in (j2-16p - (j2T9p . We test the validity of the eikonal and 
adiabatic approximations used in the derivation of F^-^^ . 
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Fig. 2. Zij integrated TEL potential as a function of bij for (a) -Enn ~ 300 MeV and (b) -Enn = 
800 MeV. The solid (dotted) curve represents the real (imaginary) part, while the dashed curve 
corresponds to the Zij integration of v{rij). 



We assume y 3> 1 and consider the accuracy of 

r^f{b,,) « 1 - exp[-i?°7(nt;rei)], 
taking t'°'^ as t'°'^. From (j2-16p through (|2-19p it is clear that 



(3-5) 



(3-6) 



is the Glauber approximation of the transition matrix of a fictitious PA collision in 
which only nucleons i in P and j in A interact. The exact wave function ip describing 
the same process satisfies the Schrodinger equation 



{K + hp + hA + r°^{r,j) -E)ij = 0, 



(3-7) 



where K = —{h\/iif/{2^a) with //q, the reduced mass of PA system, rij = R+Xi—yj 
and E = E-px + ep + ca with the incident energy E-px and the intrinsic energy ep 
(ca) of P (A). The exact transition matrix T^x calculated using ip is 



loc/ 



where 



Tex = t'°7r,,) + t^°7r,,)GoTex, 



Gn 



E-K 



(3-8) 



(3-9) 



/ip — /lA + ie 

First, we test the validity of the eikonal approximation. For this purpose, we 
apply the adiabatic approximation to Eq. ()3-7p . The equation is reduced to 



{K + t 



loc/ .^ 
ij V «J; 



EvA)i^ 



AD 







(3-10) 
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where the adiabatic approximation h-p + /iA has been replaced by the ground state 
energies ep + eA- We then obtain the transition matrix Tad under the adiabatic 
approximation as 

(3-11) 



^AD ^ tl°c(^^^.) +iloc(^^^.)cAD^AD^ 



where 



^0 



AD 



1 



^PA -K + ie 



(3-12) 



Since xi and yj are simply parameters in p-lOp . we can regard K as K = — (^Vr^^ )^/(2/io 
and then solve (j3-10p without the eikonal approximation. 

Now it is possible to test the accuracy of the eikonal approximation by comparing 
Tq\ with Tad- As an example, we consider ^He+^'^^Pb scattering. As mentioned in 
§ 13.11 t^°'^{rij) is dependent on 9. In order to avoid this complication, we use the 
^-independent (spherical) potential 



*0 Vij) 



1 1 d 

71 Tij QjTij 



cxD rioc 



t'°%hAh,dh 



ijj^ij^'^ij 



(h: 



J.2 \l/2 



(3-13) 



instead of t {vy ) . As found from the discussion in § 12.31 the two potentials tg (''"ij ) 
and t {vij) give the same transition matrix under the eikonal approximation. 

Figure [3] shows the scattering amplitudes, /gi(q) = —{2'ir)'^fi/h? x Tqi and 
/AD(q) = -(27r)^/x//i^ X Tad, of a single NN collision in ^He+^°^Pb scattering at 
the laboratory energy of E'pa = 300A^p MeV. The solid (dashed) and dotted (dash- 
dotted) curves show, respectively, the real and imaginary parts of /ad(q) (/gi(q))- 
The agreement between calculations with and without the eikonal approximation is 
excellent. 




Non-eikonal (real) 
Non-cikonal (imaginary) 
Eilconal (real) 



— ■ — ■ Eilconal (imaginary) 



4 8 12 16 20 24 28 32 

q (fm-') 

Fig. 3. On-shell scattering amplitudes, /gi('7) and fAoiQ), of a single NN collision in ^He+^^^Pb 
scattering t]^'^{rij) at the laboratory energy of EpA ~ SOOA^'p MeV. 
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The denominator of Gg in the momentum representation is 

^ -hvrjq,--^ , (3-14) 



where t^rei = hk/fia, EpA = h'^k /2iia, q = k — k' and q^ is the z component of q. 
In the eikonal approximation, the q^ term in the denominator is neglected. This is 
reahzed in the large-^^ hmit with v^^i fixed. In general, PA scattering has a larger 
fia than NA scattering, and //q, becomes minimum for NN scattering. Thus, the 
eikonal approximation is generally better for PA scattering than for NA scattering 
but becomes worse for NN scattering, when the scatterings of different systems with 
common v^ei are compared with each other. 

In order to confirm the validity of the adiabatic approximation, we compare 
the solution T^x of Eq. (j3-8p with T^^ of Eq. (j3Tip for nucleus-proton scattering 
in which the nucleus is described by a core+nucleon (c + n) two-body system. We 
consider the two reactions ^^Be+p and ^^Ca+p at an incident energy of 300 MeV 
per nucleon. We use a Woods-Saxon interaction between c (^'^Be or ^^Ca) and n 



Kn(r) = l^o(l + exp^-^j , (3-15) 

where r is the displacement of n from c. The potential parameters and the resulting 
neutron separation energy Sn are shown in Table [H We solve the three-body scatter- 

Table I. Potential parameters in Eq. (|3-15p and neutron separation energy Sn for ^°Be-n and ^®Ca- 





Vo (MeV) 


Ro (fm) 


ao (fm) 


S„ (MeV) 


'°Be-n 


51.60 


2.996 


0.52 


0.503 


^^Ca-n 


51.68 


4.343 


0.67 


15.64 



ing problem by the method of continuum-discretized coupled channels (CDCC)^^''^^^ 
with and without adiabatic approximation to the c-n internal Hamiltonian, i.e., hp 
in Eq. (|3-7p . For simplicity, we neglect the internal degrees of freedom of the target 
proton; we then have /ia = in the present case. The maximum r is 60 fm, and s-, 
p-, d- and /-waves for the c-n relative wave function are included. The momentum 
bin is truncated at 1.5 fm~^ for each partial wave and divided into 30 (15) discretized 
states for the s-wave (p-, d- and /-waves). The coupled-channel potentials are calcu- 
lated using to°'^(?'jj) of Eq. (j3T3p . Note that the interaction between c and the target 
proton is not included since we solve the Schrodinger equations (|3-7p and (|3-10p . The 
relative wave functions ifj and ip are integrated up to 20 fm and connected to the 
standard asymptotic form. 

In Fig. m we show the total breakup cross sections for ^^Be+p (left panel) and 
'^'^Ca+p (right panel) at 300 MeV per nucleon as a function of the center-of-mass 
scattering angle of the c-n two-body system. In each panel the solid and dashed lines 
indicate the results with and without adiabatic approximation, respectively. One sees 
clearly that the adiabatic approximation works very well for not only ^^Be with small 
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Fig. 4. Total breakup cross sections for ^^Be+p (left panel) and ^"Ca+p (right panel) at 300 MeV 
per nucleon. The horizontal axis is the scattering angle of the center of mass of the fragments 
c and n. The solid and dashed lines show the results of CDCC calculation with and without 
adiabatic approximation, respectively. 



Sn but also ^'^Ca with quite large Sn- Thus, in the energy region considered here, 
namely, a few hundred MeV per nucleon, the adiabatic approximation and hence 
Eq. p-lOp is shown to be valid. 

3.3. Relationship between present theory and conventional Glauber calculations 

In this subsection, we clarify the relationship between the present theory and the 
conventional Glauber calculations using the empirical NN profile function (jl-4p '^) or 
using the modified profile functions including the in-medium effects.^) For simplicity, 
we consider the case of Y ^ 1, and the superscript loc of t^°^ is suppressed. 

In the present theory, the transition matrix of the PA scattering in coordinate 
representation is 



T = C 



1 - JJexp(-zfjj//i 



Vrel, 



«J 



The corresponding matrix in the Glauber theory is 



Tgi = C 



n(i 



iifl'/hv.el) 



(3-16) 



(3-17) 



«j 



where tfj" 



— i/lVrei-^NN using F^^ of (|1-4|) . Since t^™ is adjusted to fit the data 
on NN scattering, it is essentially identical to tij. One sees, therefore, that the 
conventional Glauber calculation with (jl-4p is justified provided that (a) r ~ t is a 
good approximation and (b) fij/hv^ei is small enough to warrant exp{—ifij/hvre\) ~ 
1 — ifij/hv-cei- At very high energies, where the replacement of r by t is valid and 
the amount of multiple scatterings by t is negligible, that is, in the limit that the 
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impulse approximation is good, the present theory agrees using the conventional 
Glauber theory with the empirical profile function F^ of (jl-4p . Conversely, if either 
of conditions (a) and (b) is not satisfied, the conventional Glauber procedure is not 
justified by the present approach. 

In the heuristic approach in Ref. 5), t|j° is replaced by the z integration of gij, 
Qij , in the conventional Glauber theory. The corresponding transition matrix of the 
PA scattering is 

r(s) = C[l - n(l - ^~9^J/hVrel)\ ■ (3-18) 

One can find from (j3-16p and (j3-18p that even if fij ~ gij, T^^) agrees with T in the 
lowest (first) order of g/hv^d, but not in higher orders. Thus, the heuristic approach 
in Ref. 5) gives a better description of PA scattering than the Glauber theory, but 
it is still not perfect, since it contains no higher order terms of g/hv^f.\. 

For the calculation of optical potential, it is common to take the optical limit in 
which only the first order of g/hv^-ei is taken into account. In this limit, the present 
theory with g as t reduces to the folding model with g, which has been successfully 
used in reproducing NA scattering data.^^"^^^ 

§4. Summary 

We present in this paper a new version of the Glauber theory for nucleus-nucleus 
collisions based on the multiple scattering theory of Watson (MST),^^^ making use of 
the KMT formalism. ^^' The input of the theory is the effective nucleon-nucleon (NN) 
interaction r of MST, which has a weaker short-range repulsion than the bare NN 
nuclear force potential, which makes the eikonal approximation much more reliable. 
We tested the validity of the eikonal and adiabatic approximations to NN collision 
using T in nucleus-nucleus scattering at the laboratory energy of £'pa = 300A^p MeV 
and showed that these approximations are good at the intermediate energies. 

At very high energies where the replacement of r by t is valid and multiple 
scattering by t is negligible, that is, in the limit that the impulse approximation is 
good, the present theory gives the same T matrix of PA scattering as the conventional 
Glauber theory using the empirical profile function F^^ of p-4p . 

When the g matrix is used as r, the present theory can take into account the 
nuclear-medium effects in the analyses of NA and PA scattering. Therefore, the 
present theory is also applicable in the intermediate-energy region. The present 
theory also provides a theoretical foundation to heuristic Glauber calculations^' in 
which the profile function is modified to reproduce the in-medium NN cross section 
calculated from the g matrix, although the heuristic calculations include no multi- 
ple scattering by g. For the calculation of optical potential, the present formalism 
reduces in the optical limit to the folding model using g matrix, which has been 
successfully used in reproducing nucleon- nucleus scattering data.^^^'^^^ 

Thus, the present theory unifies standard methods such as the conventional 
Glauber method^-* using the empirical NN profile function (jl-4p .'^-* heuristic Glauber 
calculations with the modified profile functions including the in-medium NN effects, ^^ 
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and the folding potential method using the g matrix.^''' ^^' In our forthcoming paper, 
we propose a way of deriving r with no local density approximation, and compare 
the result with those obtained using the t and g matrices. 

Acknowledgements 

This work has been supported in part by a Grant-in- Aid for Scientific Research 
(18540280) from the Ministry of Education, Culture, Sports, Science and Technology 
of Japan. 

References 

1) I. Tanihata et al, Phys. Rev. Lett. 55 (1985), 2676; Phys. Lett. B 287 (1992), 307. 

2) R. J. Glauber, Lectures in Theoretical Physics (Interscience, New York, 1959), Vol. 1, p. 
315. 

3) J. S. Al-Khalili and J. A. Tostevin, Phys. Rev. Lett. 76 (1996), 3903. 

J. A. Tostevin and B. A. Brown, Phys. Rev. C 74 (2006), 064604, and references cited 
therein. 

4) Y. Ogawa, T. Kido, K. Yabana and Y. Suzuki, Prog. Thcor. Phys. Suppl. No. 142 (2001), 
157, and references cited therein. 

5) C. Xiangzhou, F. Jun, S. Wenqing, M. Yugang, W. Jiansong and Y. Wei, Phys. Rev. C 
58 (1998), 572. 

A. de Vismes, P. Roussel-Chomaz and F. Carstoiu, Phys. Rev. C 62 (2000), 064612. 
R. E. Warner, L J. Thompson and J. A. Tostevin, Phys. Rev. C 65 (2002), 044617. 

6) W. Horiuchi, Y. Suzuki, B. Abu-Ibrahim and A. Kohama, Phys. Rev. C 75 (2007), 044607. 

7) R. J. Glauber and G. Matthiae, Nucl. Phys. B 21 (1970), 135. 

8) R. B. Wiringa, V. G. J. Stoks and R. Schiavilla, Phys. Rev. C 51 (1995), 38. 

9) S. J. Wallace, Ann. of Phys. 78 (1973), 190; Phys. Rev. D 8 (1973), 1934, and references 
therein. 

10) C. W. Wong and S. K. Young, Phys. Rev. C 12 (1975), 1301. 

S. K. Young and C. W. Wong, Phys. Rev. C 15 (1977), 2146; Phys. Lett. B 77 (1978), 41. 

11) S. J. Wallace, Adv. Nucl. Phys. 12 (1981), 135, and references cited therein. 

12) K. M. Watson, Phys. Rev. 89 (1953), 575. 

13) A. K. Kerman, H. McManus and R. M. Thaler, Ann. of Phys. 8 (1959), 551. 

14) N. M. Queen, Nucl. Phys. 55 (1964), 177. 

D. R. Harrington, Phys. Rev. 184 (1969), 1745. 

J. M. Eisenberg, Ann. of Phys. 71 (1972), 542. 

V. B. Mandelzweig and S. J. Wallace, Phys. Rev. C 25 (1982), 61. 

15) G. Takeda and K. M. Watson, Phys. Rev. 97 (1955), 1336. 

16) A. Picklesimer and R. M. Thaler, Phys. Rev. C 23 (1981), 42. 

17) W. G. Love and M. A. Franey, Phys. Rev. C 24 (1981), 1073 
M. A. Franey and W. G. Love, Phys. Rev. C 31 (1985), 488. 

18) L. Ray, Phys. Rev. C 41 (1990), 2816. 

L. Ray, G. W. Hoffmann and W. R. Coker, Phys. Rep. 212 (1992), 223, and references 
cited therein. 

19) H. F. Arellano, F. A. Brieva and W. G. Love, Phys. Rev. C 52 (1995), 301. 

20) G. R. Satchler, Direct Nuclear Reactions (Oxford University Press, New York, 1983), and 
references cited therein. 

21) G. Bertsch, J. Borysowicz, H. McManus and W. G. Love, Nucl. Phys. A 284 (1977), 399. 

22) J.-P. Jeukenne, A. Lejeune and C. Mahaux, Phys. Rev. C 16 (1977), 80; Phys. Rep. 25 
(1976), 83. 

23) F. A. Brieva and J. R. Rook, Nucl. Phys. A 291 (1977), 299; ibid. 291 (1977), 317; ibid. 
297 (1978), 206. 

24) N. Yamaguchi, S. Nagata and T. Matsuda, Prog. Theor. Phys. 70 (1983), 459. 
N. Yamaguchi, S. Nagata and J. Michiyama, Prog. Theor. Phys. 76 (1986), 1289. 

25) L. Rikus, K. Nakano and H. V. von Geramb, Nucl. Phys. A 414 (1984), 413. 



A New Glauber Theory based on Multiple Scattering Theory 17 

L. Rikus and H. V. von Geramb, Nucl. Phys. A 426 (1984), 496. 

26) K. Amos, P. J. Dortmans, H. V. von Geramb, S. Karataglidis and J. Raynal, Adv. Nucl. 
Phys. 25 (2000), 275. 

27) T. Fummoto and Y. Sakuragi, Phys. Rev. C 74 (2006), 034606. 

28) M. Kamimura, M. Yahiro, Y. Iseri, Y. Sakuragi, H. Kameyama and M. Kawai, Prog. 
Theor. Phys. SuppL No. 89 (1986). 

N. Austern, Y. Iseri, M. Kamimura, M. Kawai, G. H. Rawitscher and M. Yahiro, Phys. 
Rep. 154 (1987), 125. 

29) N. Austern et al., Phys. Rev. Lett. 63 (1989), 2649; Phys. Rev. C 53 (1996), 314. 



